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We introduce a statistical physics model for opinion dynamics on random networks where agents 
adopt the opinion held by the majority of their direct neighbors only if the fraction of these neighbors 
exceeds a certain threshold, p u - We find a transition from total final consensus to a mixed phase 
where opinions coexist amongst the agents. The relevant parameters are the relative sizes in the 
initial opinion distribution within the population and the connectivity of the underlying network. 
As the order parameter we define the asymptotic state of opinions. In the phase diagram we find 
regions of total consensus and a mixed phase. As the 'laggard parameter' p u increases the regions of 
consensus shrink. In addition we introduce rewiring of the underlying network during the opinion 
formation process and discuss the resulting consequences in the phase diagram. 
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PACS numbers: 89.75.Fb, 87.23. Gc, 05.90.+m 



INTRODUCTION 

Many decisions of human beings are often strongly in- 
fluenced by their social surroundings, e.g. the opinion 
of friends, colleagues or the neighborhood. Only a few 
types of decisions in few individuals emerge from abso- 
lute norms and firm convictions which are independent 
of the opinion of others. Much more common is the situ- 
ation where some sort of social pressure leads individuals 
to conform to a group, and take decisions which minimize 
conflict within their nearest neighborhood. For example, 
if a large fraction of my friends votes for one party, this 
is likely to influence my opinion on whom to vote for; 
if I observe my peers realizing huge profits by invest- 
ing in some stock this might have an influence on my 
portfolio as well; and if the fraction of physicist friends 
(coauthors) publishing papers on networks exceeds a cer- 
tain threshold, I will have to reconsider and do the same; 
the social pressure would otherwise be just unbearable. 
Lately, the study of opinion formation within societies 
has become an issue of more quantitative scientific in- 
terest. In first attempts agents were considered as sites 
on a lattice, and opinion dynamics was incorporated by 
the so-called voter model (VM) [l[ (only two neighbors 
influence each other at one timestep), the majority rule 
(MR) @, [3] (each member of a group of odd size adopts 
the state of the local majority), or the Axelrod model Q 
(where two neighbors influence themselves on possibly 
more than one topic with the objective to become more 
similar in their sets of opinions) . Imposing regular lattice 
structure on social environments is convenient however, 
most observed structures of real-world networks belong 
to one of three classes: Erdds-Renyi (ER) scale- free 
Q or small- world networks ■ This has been accounted 
for the VM J" 
topologies [111 11 
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as well as for the MR on different 
. For a review of further efforts in 
15] and citations therein. Another 
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FIG. 1: Update process for two different configurations of 
neighbors and an update threshold of p u = 0.8. The node in 
the center gets updated, (a) Three out of five neighbors are 
in a different state, so the threshold is not exceeded and the 
node stays unchanged, (b) Four out of five neighbors are in a 
different state; 4/5 > 0.8 thus the node adopts the state. 



approach to model social interaction was developed out 
of the notion of catalytic sets [r|, leading to an una- 
nimity rule (UR) model 17] on arbitrary networks in an 
irreversible formulation. 

As a realistic model for many real world situations here 
we present a reversible generalization to the UR and MR 
models introducing an arbitrary threshold governing up- 
dates ('laggard' parameter). The UR and MR are ex- 
tremal cases of the model. In [r| the idea of a threshold 
was introduced in the context of investigating the origin 
of global cascades in ER networks of 'early-adopters'. In 
contrast to this work, where updates were only allowed 
in one direction, i.e. irreversible, the following model is 
fully reversible in the sense that two opinions compete 
against each other in a fully symmetric way. 
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THE MODEL 

Each individual i is represented as a node in a net- 
work. The state of the node represents its opinion on 
some subject. For simplicity we restrict ourselves on bi- 
nary opinions, yes/no, 0/1, Bush/Mother Theresa, etc. 
Linked nodes are in contact with each other, i.e. they 
'see' or know each others opinion. The opinion forma- 
tion process of node i is a three-step process (see FigH]): 
Suppose i is initially in state 'O'('l'). 

• Check the state of all nodes connected to i. 

• If the fraction of state 'l'('0')-nodes of i's neighbors 
exceeds a threshold p u , i adopts opinion 'l'('O'). 

• Otherwise i remains in state 'O'('l'). 

As a substrate network we chose random graphs [5j, 
i.e. N nodes are randomly linked with L links (self- 
interactions are forbidden), the average connectivity be- 
ing k = L/N. The update threshold necessary for a 
node's change of opinion, p u has to be higher than 0.5 in 
order to be meaningful in the above sense. The update is 
carried out asynchronously. In a network containing N 
nodes, at time t, there are A® nodes with opinion '0' and 
A] nodes with opinion '1'. The relative number of nodes 
are a° = A°J X /N. One time step is associated with 
applying the update procedure N times, i.e. each node 
gets updated once per timestep on average. As time goes 
to infinity, the relative population of nodes with opinion 
0/1 will be denoted by a^ 1 . 

ANALYTICAL AND NUMERICAL RESULTS 

To derive a master equation for the evolution of this 
system we calculate opinion-transition probabilities via 
combinatorial considerations in an iterative fashion, mo- 
tivated by [lq . A master equation for a° is found ex- 
plicitly, the situation for a\ is completely analogous. At 
t = 0, we have a fraction of aj] nodes in state '0'. The 
probability that at time t one node belonging to a° will 
flip its opinion to '1' is denoted hyp® . This probability 
is nothing but the sum over all combinations where more 
than a fraction of p u of the neighbors are in state '1', 
weighted by the probabilities for the neighboring nodes 
to be either from a° or a\ = (l — a°), 

PV l = t Q)(l-a?)>°)^ , (1) 

where [.] denotes the ceiling function, i.e. the near- 
est integer being greater or equal. The same consid- 
eration leads to an expression for the opposite tran- 
sition pl^°, where 1 and are exchanged in Eq.([T]). 



The probability for a node to be switched from '0' to 
'1', Aq^ 1 , is the product of the transition probability, 
Pt^ 1 , and the probability to be originally in the frac- 
tion Op, i.e. Aq^ 1 = Pq^clq. The same reasoning gives 
Aq^° = pl^° (l — dp) and provides the master equation 
for the first time step (i.e. updating each node once on 
average) , 

= + *tr° Ar 1 . (2) 

Let us now examine some special cases. 

The low connectivity limit 

For sufficiently low connectivities k there are no more 
possible updates in the network after the first iteration. 
For a given update threshold p u , this is the case if a 
change in opinion requires all neighboring nodes to have 
the same state. To see this more clearly, consider the case 
of a network with constant k — 2 (ID circle). Choose a 
node whose state is e.g. '0'. There are four possible 
configurations of neighbors: both being in state '0', one 
being in '0' and the other in '1' and both being in '1'. 
Irrespective of p u , only the latter configuration allows 
an update. For all other cases at least one neighbor in 
state '0' must be updated to '1', i.e. has to have two 
neighboring nodes in '1'. But this is not possible, since 
at least one neighbor will always be in '0'. The same 
holds for higher values of fc, as long as every neighbor 
has to hold the same opinion to allow an update, i.e. 
we effectively use an unanimity rule. For the special case 
k = 2 the final population in state '0' is given by a^, = a° . 
Inserting this in Eq.([2]) yields 

a° oc =3(aS) 2 -2(aS) 3 . (3) 

A comparison between the theoretical prediction of 
Eq.Q and the simulation of this system (on a regular 
ID circle network with N = 10 4 ) is seen in Figj2ta). 

Higher connectivities 

For higher connectivities there are much more config- 
urations allowing for potential updates, the evolution 
does not stop after one single iteration. At the sec- 
ond time step the update probability is given by the 
product of the transition probability at t = 1 and the 
probability to be initially in the respective state, re- 
duced by the probability to already have undergone this 
transition during the first time step. We have, for ex- 
ample, A}^° = (p\^° -Po^°) (l - a o)- For arbitrary 
times t this is straight forwardly seen to be A^° = 
(pt - * — Pt-x) (l — a>o)i an d the master equation be- 
comes a° +1 = a° + Aj^° — Aj - * 1 . Inserting for in 
a recursive way yields the master equation 

a t % 1 =a°+ P r , (l-a , )-*r i «S ■ (4) 
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FIG. 2: Asymptotic population sizes of the 'O'-state fraction, a^, as a function of its initial size, ajj, for N = 10 4 , p u = 0.8. (a) 
k = 2 for all nodes (ID circle), (b) ER graph with k = 9000 and (c) ER graph with k = 10. 



Again, theoretical predictions of Eq.((3]) agree perfectly 
with numerical findings, Fig[^b). Three regimes can be 
distinguished: two of them correspond to a network in 
full consensus. Between these there is a mixed phase 
where no consensus can be reached. 

High connectivity limit. For the fully connected net- 
work the asymptotic population sizes can easily be de- 
rived: if <2q > p u or <2q < 1 — p u consensus is reached. 
For 1 — p u < Oq < p u the system is frustrated and no up- 
date will take place, giving rise to a diagram like FigfJIb). 
Compared to FigJ^fa) a sharp transition between the con- 
sensus phases and the mixed phase has appeared. We 
now try to understand the origin of this transition. 

Intermediate regime. The transition between the 
smooth solution for the final populations as a function of 
a{j and the sharp one for higher connectivities becomes 
discontinuous when the possibility for an individual node 
to get updated in a later timestep ceases to play a neg- 
ligible role. Systems with small update probabilities will 
then be driven towards the consensus states. However, 
if the initial populations are too far from the consensus 
states they will not be reached. For p u = 0.8 the sharp 
transition arises for values of k around 10. FigJ^c) shows 
simulation data for ER graphs with N = 10 4 nodes and 
k = 10 with p u — 0.8. Here we already find two regimes 
with consensus and an almost linear regime in-between. 
The analytical curve obtained from numerical summa- 
tions of Eq. (j4|) resembles the qualitative behavior of the 
simulations up to finite-size deviations. The dynamics of 
the system is shown in the phase diagram, Fig[3]Ja). It 
illustrates the size of the respective regimes and their de- 
pendence on the parameters a° and connectedness k/N . 
The order parameter is a^. Along the dotted lines a 
smooth transition takes place, solid lines indicate dis- 
continuous transitions from the consensus phase to the 
mixed phase. The change from smooth to sharp appears 
at k/N w 0.01. For larger p u the regions of consensus 
shrink toward the left and right margins of the figure. 

So far we assumed static networks. However, this is 
far from being realistic, as social ties fluctuate. We thus 



allow links to get randomly rewired with the rewirement 
process taking place on a larger time scale than the opin- 
ion update, since otherwise the new connection would not 
lead to state changes. Let us assume that the number of 
rewired links per rewirement-timestep is fixed to L' , so 
that it becomes natural to define a social temperature, 
T = L'/L. T quantifies the individual's urge to recon- 
sider a topic with new acquaintances, or equivalently, the 
fluctuation of ties in their social surrounding. 

The evolution of opinions in a network at T ^ is 
as follows: We fix a network and perform the same dy- 
namics as for T — 0, until the system has converged and 
no further updates occur. Then perturb the system by 
a rewirement step and randomly rewire L' links among 
the N nodes (N and L are kept constant over time), in- 
crease the time-unit for the rewirement steps by one and 
let the system relax into a (converged) opinion config- 
uration. Iterate this procedure. Note that this process 
can be viewed as a dynamical map of the curves shown 
in Figs(21Ja)-(c). With this view it becomes intuitively 
clear that consensus will be reached for a wider range 
of parameters, where the time to arrive there crucially 
depends on the value of k. 

To incorporate the temperature effect in the master 
equation we introduce the second timescale and denote 
the population in state '0' as a°p Here t is the time for 
the update process as before and t is the time step on the 
temperature time scale, i.e. counts the number of rewire- 
ment steps. We use a,g = % . can be obtained 
from a°, )0 = lin^ (a%_+ A 1 ^ -A^ 1 ) for high k, 
and from Eq.([2]) for low k, when we only observe up- 
dates during the first iteration. This evolution is nothing 
but a dynamical map. The probabilities to find a con- 
figuration of neighbors allowing an update are no longer 
given only by A^ 1 and A^ , instead we have to count 
the ones constituted by a rewiring, which happens with 
probability T. That is why we can consider this kind 
of evolution as a dynamical map of the former process, 
with o as the initial population for the first rewire- 
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FIG. 3: Phase diagram for a function of initial frac- 

tion size a® and connectedness, k/N. Simulations where per- 
formed with ER graphs with N — 10 3 and p u = 0.8. Two 
symmetrical regions of consensus and a mixed phase in be- 
tween are observed. The dotted line indicates a smooth tran- 
sition, the solid line a discontinuous one. Inset: Detail for 
small k. Arrows mark the change from smooth to sharp tran- 
sitions, positioned at k/N ~ 0.01. (b) Phase diagram for 
a tL,oo- Technically adjacency matrices with N = 10 4 were 
generated and checked by Monte-Carlo simulations whether 
they allow an update at fixed a[] and k. 



merit step evolving to 1 , and so on. The transition 
probabilities are now given by TA^° and TA]f , since 
only new configurations can give rise to an update. We 
thus assume the master equation for a system at T ^ 
after the first rewiring to be 



"oo,t+i 



lim (a°- + T(Ah>° 



(5) 



Furthermore, one expects the existence of a critical value 
k c , below which the intermediate regime (mixed state) 
will disappear. This will occur whenever there is no 
chance that a configuration of neighbors can be found 
leading to an update. The value for k c can be easily 
estimated: Say we have a node in state T' and ask if 
an update to state '0' is possible under the given circum- 
stances. For a given k this requires that there are at least 
\kpu] neighbors in state '0' present in the set A®. If k is 
above the critical value k c it occurs that even if all nodes 
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FIG. 4: Half-life time r (i until half the population reached 
consensus) vs. relative number of neighbors for T = 
0.25,0.5, 1. Inset: Same in log-log scale. Scaling around the 
pole k c /N ~ 0.61 with an exponent 7 ~ 7.4 is suggested. 10 3 
initial populations with oq = 0.5 and N = 10 2 were averaged. 



from Aq were neighbors of the node in state T', there are 
still too many other neighboring nodes (which are then 
necessarily in state T') to exceed the update threshold. 
This means that we can not have updates if \kp u ~\ > Aq, 
and we get 
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For p u = 0.8 and = 0.5, k c « 0.617V. We next con- 
sider the time-to-convergence in the system. To this end 
we measure the half-life time r, of initial populations at 
a[] = 0.5 for different connectivities k, see Figj4] The 
figure suggests that the observed scaling of r could be of 

power-law type, with a pole at k c /N, i.e. r oc 

The estimated critical exponent 7 « 7.4 seems to be in- 
dependent of temperature. Note, that the estimate is 
taken rather far from the pole at k c , which suggests to 
interpret the actual numbers with some care. 

The phase diagram for the T / system is shown 
in Figj3](b). There are still three regimes, which are ar- 
ranged in a different manner than before. Consensus is 
found for a much wider range of order parameters; the 
mixed phase is found for high connectivities, i.e. k > k c . 
The value of k c at et[] = 0.5, as found in FigOJb), is 0.63, 
slightly above the prediction of 0.61. This mismatch is 
because we used networks with inhomogeneous degree 
distributions (Poisson) . Whether a network allows for an 
update or not is solely determined by the node with the 
lowest degree k, which explains why we can still observe 
updates when the average degree k is near to but already 
above k c . Systems in the mixed phase are frustrated. k c 
is linear in a[J which we confirm by finding a straight line 
separating the frustrated drom the consensus phase, see 
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FiglS^b). For larger p u the regions of consensus shrink. 

CONCLUSION 

Summarizing we presented a model bridging the gap 
between existing MR and UR models. Opinion dynamics 
happens on static random networks where agents adopt 
the opinion held by the majority of their direct neighbors 
only if the fraction of neighbors exceeds a pre-specified 
laggard-threshold, p u . The larger this parameter the 
more stimulus the agent needs to adapt his opinion to 
the one of his direct neighborhood. This system shows 
two phases, full consensus and a mixed phase where opin- 
ions coexist. We studied the corresponding phase di- 
agram as a function of the initial opinion distribution 
and the connectivity of the underlying networks. As the 
laggard-parameter p u increases the regions of full con- 
sensus shrink. We introduced rewiring of the underlying 
network during the opinion formation process and discuss 
the resulting consequences for the phase diagram. This 
social temperature introduced here differs from the usual 
temperature of statistical mechanics. It accounts for link 
fluctuations and not for the fluctuations of the state of 
the nodes. For T > 0, the system can escape the frozen 
state aH 1 / 1, and global consensus can be obtained. In 
the case of usual temperature (opinions of nodes switch 
randomly) a different behavior is expected. For low 
temperature, the system also can escape the frozen state, 
however for higher values of T the system undergoes a 
transition from an ordered to an unordered phase, where 
= 1/2. Even though laggards sometimes enjoy a 
bad reputation as being slow and backward-oriented, so- 



cieties of laggards are shown to have remarkable levels of 
versatility as long as they are not forced to interact too 
much. 
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